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Abstract 



^ We investigate the statistical properties of two dimensional random cellular systems (froths) in 
term of their shell structure. The froth is analyzed as a system of concentric layers of cells around 
\^ a given central cell. We derive exact analytical relations for the topological properties of the sets of 
^ cells belonging to these layers. Experimental observations of the shell structure of two-dimensional 
^ soap froth are made and compared with the results on two kinds of Voronoi constructions. It is 
^ found that there are specific differences between soap froths and purely geometrical constructions. 
In particular these systems differ in the topological charge of clusters as a function of shell number, 
in the asymptotic values of defect concentrations, and in the number of cells in a given layer. We 
derive approximate expressions with no free parameters which correctly explain these different 
behaviors. 



1 Introduction 

Materials consisting of cellular structures such as metal grains and biological tissues are common 
in nature [1, 2]. Among these systems, soap froth is considered as the paradigm for the study of 
trivalent two-dimensional cellular structures. The structural analysis of two-dimensional cellular 
patterns have been made by many researchers and many interesting results have been obtained 
[3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. These studies of the topological properties of soap froth can 
be summarized in several laws such as the Lewis' law[13] on the the statistics of cell area, the 
von Neumann's law[14] on the growth rate of n-sided cells, the scaling law[6] on the probability 
distribution of the cells, and the Aboav-Weaire's law on nearest neighbor correlation [15, 16]. So 
far, the evolution of soap froth after the scaling state defined as stationary probability distribution 
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of cell sides have been explained by several theories [17, 18, 19, 20], indicating that correlation 
effects are not manifested in the analysis of area scaling law. However, more detailed analysis 
beyond area scaling law has been done[21] and suggests strongly the importance of clarifying 
the role of correlation effects. Moreover, there has never been any experimental investigation till 
recently[22] on the statistical properties and correlation effects beyond nearest neighbors in two- 
dimensional soap froth. With the advance of experimental techniques and data analysis[23, 24, 22], 
it is natural to investigate the structural characteristic of soap froth beyond nearest neighbors. 
Here we derive systematically exact expressions for the topology of the froth structure beyond 
first neighbours. Our results will serve as a mathematical framework for our data analysis, which 
enhances our ability to test the current theoretical ideas and our understanding of topological 
ordering processes in soap forth [25, 26, 27, 28, 29]. Our aim is to extract more information 
on the differences between physical soap froths and computer generated two-dimensional cellular 
patterns, such as the Voronoi construction of points randomly scattered on the plane. 

Soap froths are trivalent (three edges meeting at a vertex) two-dimensional cellular pattern. Any 
two-dimensional {2D) froth can be analyzed as structured in concentric layers of cells around a 
given central "germ" cell[25, 30, 22]. The first layer of the shell structure consists of the neighbours 
of the germ cell, the cells of the second layer are the neighbours that are bounding externally the 
first layer, etc. (see Fig.(l)). More precisely, the layers are closed rings of cells which are at the 
same topological distance from the germ cell (the topological distance between two cells is defined 
as the minimum number of edges that a path must cross to connect these two cells). The shell 
is closed loop of edges which are the interface between two subsequent layers (see Fig.(l)). (The 
definition here is consistent with ref. [30]. Note that in our previous work[22], the word shell is 
used to indicate the layer used in this paper. ) 

The shell structure has some physical relevance in natural froths and in problems of topologically- 
stable partition of space by cells. In foams the evolution of the cellular structure is ruled by the 
process of diffusion of the gas inside a given bubble towards its neighbours. The topological dis- 
tance j between two cells is the minimum number of soap membranes that the gas molecules have 
to cross to pass from one cell to the other. The molecules which diffuse from a given central bubble 
move through the shell structure reaching the cells of each layer with approximately equal prob- 
ability. An analogous problem where the shell structure has a physical relevance is the random 
walk through a barrier network. In this case the topological distance j is the minimum number 
of barriers that must be crossed from a starting-cell to a given final cell. The number of cells in 
the layer j is the number of final states with approximately equal probabilities. More generally, 
any perturbation on a given cell (cell growth, cell division, mechanical stress, cell-coalescence, 
electrical signal...) propagates in the whole system through the shell network and reaches with 
approximately equal intensity to all cells at the same topological distance from the perturbed 
cell at about the same time. These considerations provide the motivations of introducing shell 
structure analysis for froth. 

The quantities investigated for a shell at topological distance j from a given n-sided germ cell 
are: the total number of cells in the layer (A.,), the average number of sides of the cells belonging 
to the layer {rrij), and the topological charge {Qj) of the cluster of cells delimited by the shell 
j. Here, the topological charge of a cell with n sides is defined as g = 6 — n and is defined 
as the sum over the topological charges of all the cells inside the cluster delimited by the shell j 
which separates the layer j from the layer j Note that all these quantities, Kj^nij^ and Qj 
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can depend on n. 

The experimental investigation of the statistical properties of natural and computer generated 
froths in term of the shell structure give some interesting results. For instance, we find that the 
number of cells in a layer at a topological distance j from the germ cell increases with the distance 
following a linear law Kj = C j-\-B with slope C ~ 9 (see also Ref.[22]). This is in contrast to simple 
geometrical considerations that the perimeter of the shell cluster increases with the radius with a 
slope approximately 2tt . Moreover, we find that the topological charge of the cluster bounded by 
the shell j is negative and decreases linearly with j . This behaviour is particularly surprising since 
the total topological charge of a froth is a constant finite quantity independent of the network 
itself and is related to the space curvature by the Euler formula associated with the Gauss-Bonnet 
theorem[31, 32]. Therefore, the set of cells belonging to a layer have peculiar statistical properties 
which are different from the one of the whole froth. The aim of the present paper is to study 
these peculiarities with an exact analytical approach, and to find approximate solutions with no 
adjustable parameters which can be compared with the experimental observations. 

The plan of the paper is as follow. In Section 2, we derive the statistical properties of the shell 
structure for a special class of froths called shell-structured-infiatable (SSI) [30]. These froths are 
particularly convenient since they can be constructed layer after layer in a recursive way according 
to an infiationary procedure. In Section 3, we study the correlation length in SSI Euclidean froths. 
Approximate expressions for and are also derived. In Section 4, the results obtained for SSI 
froths are generalized to non-SSI systems. In section 5, we find approximated expressions for 
and in the general case of non-SSI froths. In Section 6, experimental results are presented and 
compared with the analytical predictions. A conclusion emphasizes the main results. In Appendix 
A, a generalization of the Weaire's sum rule [16] is derived. Einally, in Appendix B, we find an 
expression for the fiuctuations of the topological charge. 

2 Statistical properties in Shell-Structured-Inflatable 
Froths 

Any froth can be analysed as structured in concentric layers of cells which are at the same 
topological distance from a given central cell. These concentric layers are the shell structure of 
the froth. In this structure the cells making the layer j can be put into two categories. Some cells 
have simultaneously neighbours in the layers j — 1 and j These cells make themselves closed 
layers and constitute the "skeleton" [sk) of the shell-structure. Other cells (or clusters of cells) are 
inclusions between the layers of the shell-skeleton (they have neighbours in the layer j — 1 or other 
topological inclusions but not in the layer j -\- 1). The shell-skeleton is itself a space-filling froth 
hierarchically organized around the germ cell. Once a germ cell is chosen, the shell-structure and 
its skeleton are univocally defined, but different germ cells generate different skeletons. We call 
a froth shell-structured-inflatable (SSI) if it is free of topological inclusions. (In this paper these 
inclusions are also called topological defects.) Eor SSI froths the shell-structure and its skeleton 
coincide. 

In this section we analyzed only SSI systems (no topological defects), the results obtained for 
such class of systems are extendible to froths with topological defects. The advantage of studying 
SSI froths is that in these systems the shell-structure can be constructed by using an inflationary 
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recursive process. In particular, one finds that the number of vertices in successive shells are 
related by the following map [30], 

where Sj = nij — 4, and nij is the average number of sides per cells in layer j. The quantity V^^ 
(Vf) is the number of vertices attached to edges directed outward from shell j (directed inward 
to shell j) (see Fig.(l)). The matrix equation (1) is the logistic map [34]. It is a dynamical map 
from the central germ cell (j = 0) to the whole froth. 

Since our 2D froth is trivalent (three edges meeting at one vertex), the number of cells in layer 
j is Kj = V~ . In term of these quantities eq.(l) can be written as a recursive equation 

Kj+i = -SjKj - Kj-i ■ (2) 

The initial conditions are Vq = Kq = and V^^ = Ki = n= number of neighbours of the central 
germ cell. 

Given the set of parameters {s^}, the solutions of the eq.(l) or (2) are particular trajectories in 
the plane [j^Kj). In the case equals to a constant s, one finds that the parameter s separates 
the map into different classes. [30] Values |s| < 2 are associated with the elliptic region which 
has bounded, finite trajectories. The region s > 2 is the hyperbolic region associated with 
exponential trajectories. The point s = 2 divides the elliptic from the hyperbolic region and 
corresponds to tilings of the Euclidean plane. In this case one finds that the solutions are linear 
trajectories oc j. In general, when s depends on j, one can note that bounded trajectories 
always correspond to froths which are tiling elliptic manifolds whereas unbounded trajectories 
which grow faster than j correspond to tilings in hyperbolic manifolds. The Euclidean space 
has trajectories between these solutions associated with spaces of opposite curvatures. They are 
unbounded trajectories which grows asymptotically as a linear law [K^ oc j). 

One finds that the total topological charge of a cluster of cells bounded by the shell j can be 
written in terms of numbers of vertices coming in and going out from the shell j 

= ^(6 - nO = 6 - V+ + V- (3) 

i 

where the sum runs over all the cells i in the cluster and rii is the number of sides of the cell 
z[35]. Equation (3) is a general expression (valid also for non SSI systems). It states that the 
total topological charge inside a cluster depends only on its boundary. Eor SSI froths eq.(3) can 
be rewritten as 

= 6-/^,+! + /^, . (4) 

The inverse of this equation provides a relation between the number of cells in the layer j and the 
topological charge of the cells inside the cluster delimited by the shell j 

3-1 3-1 

Kj = QJ-J2Q^ = QJ-Q + n-J2 5 (5) 

8=0 8=1 

where we define Qq = 6 — n. Let us now make use of these recursive relations in order to evaluate 
the quantities Kj^ Qj, rrij^ in the shell structure. 
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2.1 First shell: the Aboav-Weaire's law 



Consider a shell-structure around an n-sided germ cell. The number of cells constituting the first 
layer (j = 1) is 

Ki = n , (6) 
the topological charge inside the cluster delimited by the first shell is by definition 

= (6 - n) + (6 - rm)K\ . (7) 

Using the sum rule (46) (see appendix A), one obtains the average value 

{Qi) = -I^2 , (8) 

with /i2 = ((n— 6)^). Here the averages ((...)) are over the cell-sides distribution: ((...)) = Zln p(^)(---)5 
with p[n) being the probability of an n-sided cell in the whole froth. One can introduce the fiuctu- 
ation part of topological charge Qi as its deviation from the average value: Ti = Qi — (Qi), where 
a-priori the fiuctuation part can be any possible function of n satisfying the condition (Fi) = 0. 
In terms of the fiuctuation, the charge Qi can be written as 

= (6-n) + (6-mi)n = -/i2 + ri . (9) 

Equation (9) gives the relationship between the number of sides (n) of a given germ cell and 
the average number of sides (mi) of the cells neighbouring this n-sided cell. In literature, such a 
relation is widely studied since Aboav founds empirically the linear law: (mi — 6)n = /i2 — a[n — 6) 
[15]. One can see immediately that the Aboav-Weaire's law is obtainable from eq.(9) by imposing 
a linear form for the fiuctuation part (i.e. Fi = ei(6 — n) with ei = 1 — a). This linear dependence 
can be interpreted in terms of the screening of the central charge Qq = 6 — n hy the charges of 
the first layer. The total screening charge in the first layer is (6 — mi)n and its deviation from 
the average is 

(6 — mi)n — ((6 — mi)n) = — a(6 — n) . (10) 

We can therefore interpretate a as a screening factor: a = 1 corresponds to a total screening of 
the internal charge, whereas a = corresponds to absence of screening. 

The Aboav-Weaire's law is generally associated with the presence of topological correlations 
between nearest neighbours cells. In Appendix B we discuss in details the effects of finite range 
correlations in froths. Let's note that an arrangement of cells free of correlation (called topological 
gas [36]) has mi = 6 -|- /i2/6 [37], which leads to an Aboav coefficient a = — /i2/6 and mi is 
independent of n. 

2.2 Generic j-shell 

The average topological charge inside the cluster delimited by shell j is 

(^?.) = E((6 - mOA;-) . (11) 

8 = 1 
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Indeed, the topological charge is an additive quantity and therefore {Q^j) is the sum of the topo- 
logical charges contained in the layers i (((6 — mi)Ki)) which are inside the cluster with 'radius' 
j ^ j)- As before, we introduce the fluctuations of the topological charge in a cluster 

r. = - (g.) . (12) 

By using the generalized Weaire identity ((6 — mj)Kj) = ((6 — n)Kj) (eq.(46) in appendix A) and 
by substituting eq.(5) into (11), we can express the portion of charge contained inside the layer j 
in term of the topological-charge fluctuations 

((6 - m,)K,) = - X:((6 - n)r.) = - 'y.{{Q - n)V,) , (13) 

8 = 1 8 = 

where we used Fq = (6 — n). The total charge of the cluster inside the shell j is therefore 

(g.) = -^23 - i:(j - ^)((6 - n)r,) = - X:(j - 0((6 - n)V,) (14) 

8=1 8=0 

Note that, in the absence of fluctuations iVi = for i > 0), eq.(13) give a topological charge per 
layer equal to — /i2 and therefore a total charge which decreases linearly in j with slope —fi2- 
Using eq.(14) and the definition of Fj, one can rewrite eq.(5) in terms of the fiuctuations 

8 = 

= 6j-(6-n)-^g, 

8 = 1 

= 6j+n-6 + /i2-2^r, + ((6-n)r,) 

8=1 8=1 

= 6j-gr.+g ^^'-^^^^-^-^^ ((6-n)r.) (15) 

8=0 8=0 ^ 

These relations for Kj , Qj and m^, are exact results valid for any SSI froth. The fiuctuations 
Fj are a-priori unknown functions subjected to the constraint (Fj) = 0. Other constraints on 
these fiuctuations come from the space-filling condition. 

3 Shell statistic of Euclidean froths 

Let us consider an Euclidean 2D froth. In such a froth the number of cells per layer must 
grow linearly with the distance in the asymptotic limit. We will show that fiuctuations in the 
topological charge (Fj ^ 0) are essential for filling space with disordered (/i2 7^ 0) Euclidean SSI 
cellular systems. Consider a system with F^ = for ^ > 1 (recall Fq = 6 — n), from the last 
equality of equation (15), one obtains that the number of cells in the generic layer j increases 
quadratically with the distance Kj oc j^fi2- Such a froth is realizable only in a space with intrinsic 
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dimension D = 3 and thus it is not two-dimensional Euclidean. In natural froths and in computer 
generated cellular system the shell structure is organized in order to keep the froth Euclidean and 
experimentally one finds that the number of cells per layer increases linearly in j after the first 
few layers[22]. Such organization can be provided by the F's which must be different from zero at 
least for the first few shells. 

We have the following interesting theorem. (See Appendix C for a proof.) Eor an Euclidean SSI 
froth with /i2 7^ that obeys the Aboav-Weaire's law, cells must be correlated at least between 
the third neighbors. This theorem can be restated into two parts: (1) If an Euclidean SSI froth 
with nonzero second moment obeys Aboav-Weaire's law, the minimum u in order that (Kj) oc j 
for j > f IS 3. (2) Under the same hypothesis as (1), and if the topological correlations vanish 
after the ^th layer, then ^ > u. Therefore, (2) also means that 1^ > 3. An important point in the 
theorem is that the Aboav parameter a is a free parameter, for otherwise, we can have an even 
smaller u. (Indeed, if we restrict a to be 1, then i/ = 1, and if we restrict a to be 2, then i/ = 2, 
and if we do not put any restriction on a, then u > 3. ) Next we will explore some consequence 
of this theorem. 

Eor a space filling Euclidean froth where the correlation length is minimal = i' = 3) and the 
Aboav-Weaire's law is satisfied with Fi = (1 — a)(6 — n), eq.(53) implies (1 — a)/i2 + ((6 — n)T2) = 
—fi2- A solution for this condition is given by r2 = (a — 2) (6 — n). Since we have = i/ = 3, we 
can set Fj = for j > 3 (from eq.(53) and (49) ). With these Fj, we can work out the solution 
for the number of cells per layer, 

Ki = n , 

K2 = 12 + /i2 + (2-a)(n-6) , 

Kj = (6 + (3-a)/i2)j-(5-2a)/i2 for j > 3. (16) 

and the topological charge 

Qi = -/i2 + (1 - a)(6 - n) , 

Q2 = -(3 - a)/i2 + (a - 2)(6 - n) , 

Qj = -(3 - a)/i2 for J > 3. (17) 

The solution eq.(16) for the trajectories Kj is qualitatively in agreement with the experimental 
data. The number of cells per layer versus j follows a linear law with slope 6 -|- (3 — a)/i2 and 
intercept —(5 — 2a)/i2- Assuming the typical values a ~ 1 and /i2 — 1-5, we obtain 9 and -4.5 
respectively for the typical slope and intercept. Note that a simple geometrical approach will 
suggest a linear growth of the number of cells of the perimeter of the shell cluster with a slope 
equal about 27r, which differs from from the one we have found. 

Equation (17) predicts that in SSI minimally correlated froths the topological charge is constant 
after the second shell. This is in contradiction with the experimental data which shows topological 
charges decreasing linearly with j. To resolve this contradiction with experimental data, we can 
consider the following scenario. If we take the fact that Qj oc j for large j, then according to 
eq. 5 implies that Kj oc j^. But our experimental froth is Euclidean with Kj oc j. Thus, the 
only possible way to get out of this contradiction with experimental data without invoking the 
results on non-SSI Euclidean froth is to conclude that the assumption that ^ = u = 3 m the above 
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analysis is incorrect. Thus this conclusion, which insists on using the results of SSI Euclidean 
froth, leads one to suspect that topological correlation has a longer range: > 3. However, this 
is actually misleading as we have no good reason to insist on the assumption that the real froth is 
SSI Euclidean. Indeed, we have enough evidence that the real froth is a non-SSI Euclidean froth, 
so that the topological correlation can vanish rather early, and still Kj oc j and Qj oc j for large 
j. This point is also valid for computer generated froth. In the next section we show that the 
topological defects provide a mechanism which correctly gives the behaviour of the topological 
charge and the slope in the linear law of the trajectories Kj. 



4 Effects of the non-inflatable inclusions 

The ideas presented in the previous sections for shell-structured froths are still applicable when 
we have non-inflatable inclusions. In this section, we indicate the appropriate corrections. We 
use the symbols n, Kj^ Qj, rrij and /i2 to indicate the quantities associated with the global froth 
(shell-skeleton plus topological inclusions) and we use the notation K'-^ , Qj^ and m**^ for the 
quantities associated with the shell skeleton only. Einally, we use K'j' and to indicate the 
number of defective cells and their average number of sides. Note that Kq = always (since any 
cell of the froth can be the germ cell). One has the simple relation 

Kj = Kf + . (18) 

In the general case when non-inflatable inclusions are present, the relations previously obtained 
for the SSI froth are applicable to the quantities associated with the shell-skeleton. In particular 
relations (4) and (5) become 

Qf = 6 - Kfl, + Kf , (19) 

and 

Kf = &j-Y.Qt . (20) 

8 = 

On the other hand, expression (3) for the total topological charge of the cluster delimited by shell 
j remains unchanged 

= 6 - y+ + V- , (21) 

here the quantities are associated to the global froth (shell skeleton plus topological in- 

clusions). This equation is the topological analog of the Gauss theorem of electrostatic. The 
charge (topological) inside a region of space is associated to the net flux (of edges) which cross 
the external surface(shell). In the absence of topological inclusion (SSI case) the flux of edges 
between two adjacent shells is uninterrupted (any edge outgoing from shell j ends in shell j -|- 1, i.e. 
V~i^i = Vj^). In the general case, the topological inclusions trap the edges-flux and the previous 
identity must be modified as follows 

y;^i = y^-v,+iKf^i , (22) 

where rj^^i is the average number of edges which are trapped by one defect in the layer j + 1. One 
can easily verify that V~ = Kj^ ■ By substituting eq.(22) and V~ into eq.(21) one obtains 

Q, = 6-K;l + Kf-r^,+,Kf^, , (23) 
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which is the generahzation of eq.(4). By substituting eq.(19) into (23) one gets 

Q, = Qf-V,+iKf+, . (24) 

Therefore, the topological charge inside shell j is equal to the charge associated with the shell- 
skeleton minus a contribution due to the defects attached to the external shell. This implies the 
important fact that the defects inside the cluster do not contribute to the total charge. Moreover, 
note that the defects attached to the external shell always decrease the topological charge in the 
shell respect to the value associated with the shell-skeleton. 

4.1 First shell and Aboav-Weaire's law 

The number of cells making the first layer around an n-sided central cell is Ki = n. By using the 
sum rule (46) obtained in appendix A one can calculate the average topological charge inside the 
first shell 

(Qi) = ((6 - n)) + ((6 - m,)K\) = ((6 - n)n) = -/i2 • (25) 

From this relation and eq.(24) follows that {Q'f') = — /i2 + (?72-^^"2)- 

As in the SSI case, one can introduce the fiuctuations of the topological charge around its 
average value 

^j = Qj-{Qj) ■ (26) 

In term of these fiuctuations and using eq.(25), the charge contained in the first layer can be 
written as 

= (6 - n) + (6 - mi)K\ = -/i2 + Ti . (27) 

This equation gives a relation between the number of sides (n) of a given cell and the average 
number of sides (mi) of the cells surrounding this n sided cell. In the particular case that 
the fiuctuation is linear in n, Fi = (1 — a)(6 — n), we arrive at the Aboav-Weaire's law (i.e. 
nnii = (6 — a)n-|-/i2 + 6a [15]). The Aboav's coefficient a can be interpreted (see §2.1) as the factor 
that represents the screening of the central charge due to the surrounding cells. By following this 
interpretation one can associate a = 1 to a total screening, and a = to the absence of screening. 
A typical value for this coefficient in natural cellular structures (soap froth, alumina cuts, etc.) 
is a ~ 1.2, values around 0.6 are characteristic of Voronoi froths constructed from Poissonian 
points, whereas random network generated by performing Tl transformations on regular lattice 
have negative Aboav's coefficients a ~ — 1. 

4.2 Generic j-shell 

The average value of the topological charge inside the layer j can be calculated by using the sum 
rule (46) (i.e. ((6 — mj)Kj) = ((6 — n)Kj), see Appendix A) and the expression (20) for Qj^, to 
obtain 

{{Q-m,)K,) = {{Q-n)Kf) + {{Q-n)Kf) 

= ((6-n)/f/)-X:((6-n)gf) 

8 = 
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J-1 



((6 - n)Kf) - Y: (((6 - n)Q,) + ((6 - nH+^Kt,)) 

8 = 

((6-n)/f/)-X:(((6-n)r.) + ((6-nK-+iAti)) • (28) 



8 = 

The average total charge of the cluster inside shell j is the sum over the charge of individual 
layers. From eq.(28) follows 

(Q^) = E (((6 - n)R%) - {j - z)(((6 - n)r.) + ((6 - n)r^,+,R%))^ . (29) 

In this expression, the i = term is ((6 — n)Kf) — j(/i2 + ((6 — n)riiKf)), since Fq = Q — n. From 
eq.(20), using eq.(29) and the definition of F^, the number of cells of the shell-skeleton in the layer 
J is 

Kf = Gj-Y.Qf 

8 = 

= 6j-'^(^F, + (g,) + r7,+iAti) 

= 6j - E (r. + ry.+iAti + (j - ^)((6 - n)Kf)^ 

+ E ~ ' ~ (((6 - ^)r8) + ((6 - nH+^Kt,)^ 

(30) 

This is the generalization of eq.(15) which takes into account non-infiatable inclusions. These 
relations have the beauty of being exact and the privilege of being useless for predicting properties 
of real space-filling cellular systems. In order to compare with experiments, we must use some 
simple physical approximations that will provide predictions for the asymptotic behaviours of A"j, 
Qj, and the percentage of defects. 



5 Euclidean froths and non SSI inclusions 

In the absence of defects and fiuctuations eq.(30) implies that the number of cells in each layer 
grows very fast [Kj oc j^). In principle, such a fast growth is realizable in the 2D plane by 
increasing continuously the roughness of the cluster-surface enlarging the available space for new 
cells in the layer. In practice, after some layers the rough surface of the cluster starts to self- 
intersect, thereby generating non SSI inclusions (topological defects in the shell structure). These 
inclusions of defective cells provide a way to smooth the shell surface and to keep the froth 
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Euclidean as a consequence. Therefore, the defects play a very important role in the froth- 
organization. 

Soap froths and computer generated random cellular systems show asymptotically a linear 
increment in the number of cells per layer with the distance from the germ cell. Let's therefore 
consider an Euclidean froth where Kj oc j for any j > u. In contrast to the SSI case, when 
we have defects, this condition of linear growth of Kj does not introduce any constraint on the 
fluctuations Fj. Indeed, the number of defects K'^ is a free parameter that the system can adjust 
in order to control the fluctuations in the topological charge and simultaneously keep the froth 
Euclidean. Analogously, the defects can play an important role in the reduction of topological 
correlations. In Section 3 we found that an SSI Euclidean froth which satisfies the Aboav-Weaire's 
law must be correlated at least between third neighbours. The defects enlarge the freedom for 
the construction of the cellular system around the germ cell by removing the constraints on 
the fiuctuations. (Without these constraints on the fiuctuations, one can construct an Euclidean 
froth which satisfy the Aboav-Weaire's law and which is correlated only between first neighbours.) 
Therefore the presence of defects can strongly reduce the correlations. 

Consider a froth which is minimally correlated and compatible with an Aboav-Weaire's law 
with a free a parameter. Erom eq.(49) (see Appendix B) it follows that such a froth can be 
uncorrelated after the first neighbours > 1). (Indeed, for the topological gas, = and this 
corresponds to a = — /i2/6 ) In this case, one gets using eq.(30) and Fi = (1 — a)(6 — n), 

K2 = 12 + /i2 + (2 - a){n - 6) - f],Kt + (1 - f]2)K^ + {{n - 6)(1 - f],)Kt) . (31) 

Transcurating the contribution from the defects in the first layer (Kf <^ 1) and fixing ri2 = 1, 
eq.(31) become 

K2 ~ 12 + /i2 + (2 - a){n - 6) . (32) 



Equation (49) gives 



By using eq.(30) one has 



(2-a)V2x 

l-a + — — j(6-n) . (33) 

IZ -\- fl2 



K3 ~ 18 + (4 - a)/i2 + (3 - 2a + %--^^) (n - 6) - (34) 

IZ -|- /i2 

where we used the same approximations as for eq.(32) plus the hypothesis rjs = 1. 

These relations are derived using the assumptions that Kf <^ 1, ^ = 2 and ri2 = rj^ = 1. The 
first condition Kf ^ 1 is quite reasonable as experimentally in real or computer generated froth, 
the first few shells have small amount of defects. The second condition is a working condition 
on the conventional wisdom that froth has rather short range correlations. In any case, this 
is a working hypothesis for deriving eq.32, 33 and 34, which are all the equations needed for 
comparison with experimental data in the next section. As for the third condition ?72 = '/s = 1, 
this is actually a reasonable assumption as rj^^i is the average number of edges which are trapped 
by one defect in the layer j -|- 1, and many defects contribute in general only one extra edge to 
the correction in the defining eq.22. 
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5.1 Asymptotic behaviour in Euclidean froths 

In an Euclidean 2D froth, the space-filling condition implies that in the asymptotic limit the 
number of cells per layer grows linearly with the topological distance and consequently at each layer 
this number is incremented by a constant amount. In a froth where the cells are uncorrelated after 
a given topological distance ^ this rate of increment must be a constant parameter characteristic 
of the whole froth and independent of the particular central cell. The dependence of Kj on the 
number of sides (n) of the central cell can only be an additive constant (i.e. Kj = Cj + B[n)). 
From a geometrical point of view this additive quantity is associated to the size of the central 
core made by the first shell and its neighbours. In a system that is uncorrelated after the first 
shell, the length of the perimeter at shell j is given by the perimeter of this core [B[n)) plus a 
term linearly dependent by the distance (Cj). The generalized Aboav-Weaire's relation (eq.(46) 
in Appendix A), gives 

((6 - mj)Kj) = ((6 - n)B{n)) for j > ^. (35) 

By assuming = 2 and using eq.(32) we obtain 

((6 - mj)Kj) = ((6 - n)K2) ^ -(2 - a)/i2 for j > 2. (36) 

Equation (36) states that the topological charge contained inside any layer is a constant quantity 
equal to —(2 — a)/i2- It follows that the total topological charge inside shell j decreases linearly 
with j 

(Qj) = —{2 — a)ii2j + const. for j ^ 1. (37) 

Note that the total topological charge in the froth must be lower than 12[31]. Indeed, a froth 
with charge Q = 12 is a closed cellular system which is tiling a surface topologically equivalent to 
a sphere. Therefore, from eq.(37) it follows that the Aboav coefficients must be smaller or equal 
than 2. 

In an Euclidean 2D froth, the number of cells per layer grows linearly with slope C asymptoti- 
cally. Let us suppose that the percentage of defects with respect to the total number Kj of cells in 
the layers is a constant A independent of the topological distance in this limit. Then the number 
of cells in the shell-skeleton must also grow linearly with j (since {K'-^) = (1 — A){Kj) and (Kj) 
is linear in j). Equation (20) indicates that {K'-^) can be linear if the average topological charge 
associated with the shell skeleton (Qj^) is constant. This implies (eq.(24) ) 

((6-m,)/^j) = (Qj-Qj-i) 

= {Qf-Qf_,)-{ri,+rKUr-mK-) 

= {{r^,Kf)-{r^,^^Kf^,)) for j » 1. (38) 
By supposing the parameter rjj = rj independent of j for large j, using eq.(38) we get 

((6 - m,)K,) ~ -vA{{K,+i) - (K,)) = -V^C for j » 1, (39) 

which predicts that the charge decreases linearly with j with a decrement of rjAC per layer. Erom 
eq.(39) and using the expression for nij obtained in Appendix B (eq.(47)) valid for layers of cells 
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uncorrelated with the germ cell (in present case for j > 2), we obtain 



g_u L^r^Qj ~ 6 _L J_ . (40) 



Supposing that (A.,) grows linearly from the second shell (i.e. fixing u = 1), from eq.(31) and 
(34) we obtain the slope 

C ~ (Ks) - {K2) ~ 6 + (3 - a)ii2 - (At) • (41) 

By comparison between eq.(36) and (39) and using expression (41) we obtain an expression for 
the percentage of defects 



(2 — a)fi2 ^ (2 — a)fi2 

6 + (3 - a)/i2 - (At) ~ 6 + (3 - a)/i2 



^ ^ ^ ^ ^ "-^^l , (42) 



where we assumed (A^) ^ 6 + (3 — a)/i2- From eq.(42) one notes that a froth can be free of defects 
only if /i2 = or a = 2. The first case (/i2 = 0) corresponds to the hexagonal lattice which is SSI 
and therefore free of defects. The second case (a = 2) corresponds to a froth where the Aboav 
parameter takes the maximum allowed value for an Euclidean froth. This is a very peculiar froth 
which has constant topological charge in the shell clusters and does not need defects to fill the 
plane. So far such a froth has not been observed. Note that a = 2 is a critical value since froths 
with a > 2 are closed elliptic systems. That might suggest that the value a = 2 is an upper limit 
which cannot be reached for Euclidean froths. 

In an earlier paper [22], we have used the shell model to test a generalization of the Aboav- 
Weaire's law to shells beyond the first. These analysis reveal a universal topological relation on 
the average number nij of sides per cell to the number of cells Kj in the j-th layer of a given center 
cell with n sides. A plot of m^Kj vs Kj shows a slope of (6 — a) for j = 1( Aboav- Weaire's law) 
and a slope of 6 for j > 2 for all samples. The results are universal for soap froths in the scaling 
state with different preparations, different times, and different temperatures. With the small rjA 
approximation and eq.40, we can now provide a quantitative explanation to the experimental 
results. The average number of neighbors nij in layer of cells uncorrelated with the central cell 
(i ^ given by eq.(47). Multiplying this expression by Kj we have 

m,K, = rr^-K, = (^6 - ^-^^j^^) (j > 0, (43) 

In the asymptotic limit, using the approximations introduced above, we obtain 

m,A^, ~ 6A^, + (2 - a)/i2 (j > 0, (44) 
which is the extension of the Aboav- Weaire's law to higher shell numbers. 



6 Experimental results and comparison with the theoret- 
ical predictions 

The soap froth chamber consists of two 1.6 cm thick retangular plexiglass plates separated by a 
0.16cm thick spacer. The effective working area for the froth is 26.7cmx36.8cm which can be filled 
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with more than twenty thousand bubbles as the starting condition. Soap bubbles of different sizes 
are pumped into the chamber through an inlet to create a relatively random froth. The chamber is 
filled with excess bubbles such that excess fiuid can be forced out of the chamber through another 
outlet after the soap froth has been drained for a few minutes by setting the chamber vertically. 
Thus the froth has a negligible volume fraction of liquid to air. The chamber is then sealed by 
closing the inlet and outlet and placed horizontally. A dark field method is used for viewing 
from above and a high resolution digital CCD camera of 1037x1344 pixels is used to capture the 
images at different stages of the evolution of the froth. The experiment starts with about 20,000 
bubbles and is about 10,000 in the scaling state defined by the stationary distribution of sides 
after about 4 hours. However, the topological charge becomes stationary at a much later time 
with about 4,000 bubbles. Runs with similar initial condition have been recorded and only data 
with stationary topological charge will be reported in this paper. 

In order to compare the physical froth with computer generated ones, we have selected two 
examples of purely geometric constructions. The first one is the Voronoi construction based on 
random Poissonian points. The second one is based on the Voronoi construction generated by 
introducing small perturbations to the triangular lattice so that we mimick the Tl transformation 
in real froth. We label the physical froth by (S), the random Voronoi froth by (V), and the 
perturbed one by (P). These two geometric constructions are two simple examples of the random 
froth and slightly ordered froth. The number of cells for these systems are: 3206 for (S), 3783 
for (V) and 9634 for (P). The values of 112 are : 1.33 (S), 1.76 (V), 1.51 (P). They all obey 
approximately the Aboav-Weaire's law with coefficient a equal respectively to 1.27ib0.05 (S), 0.69 
± 0.03 (V), 0.95 ± 0.07 (P). 

In Fig. (2) the number of cells per layer Kj is shown as function of the distance j for soap, 
Voronoi construction on set of random points in the plane and on perturbed triangular lattice. 
For all systems this number increases linearly with the distance with slopes respectively equal to 
9.45±0.1 (S), 11.0±0.2 (V), and 9.91±0.08 (P). These linear behaviours indicate that these froth 
are Fuclidean. Using the expressions for the minimally correlated Fuclidean froths obtained in 
the previous paragraph, one can predict the slopes 8.3ib0.3(S), 10.1ib0.4(V) and 9.1ib0.7(P). 

The ratios A of the defective cells respect to the total number of cells in a layer are very small 
in the second layer (0.0173(S), 0.05(V) and 0.05(P)), then increase until about the 15*'^ layer, 
to stabilize finally at an asymptotic value equal to O.lOibO.Ol (S), 0.15ib0.01 (V), and 0.13ib0.01 
(P) respectively. Using the expression (42) one can predict for rjA the values 0.136 (S), 0.23 (V) 
and 0.174 (P) respectively. We measured independently the value of rj which gives 1.1 it 0.2 (S), 
1.3ib0.1 (V), and 1.2ib0.1 (P) respectively. We deduce that the value of defect concentration A 
using these values of r] and Fq.(42) are A = 0.11 ± 0.03 (S), 0.17±0.02 (V), and 0.15±0.03 (P) 
respectively. This compares reasonably well with the measured values of A. 

Fig. (4) shows the total average charge (Q) / f^2 v.s. j. This quantity decreases with the distance, 
with an almost linear behaviour for j > 2. The measured slopes are respectively -0.7ib0.1 (S), 
-1.6±0.1 (V), and -1.35±0.05 (P). Using equation (37) one predicts -0.73±0.03 (S), -1.31±0.05 
(V), and -1.05±0.07 (P). 

The topological charge contained inside the cluster delimited by the generic shell j is an impor- 
tant physical parameter. In the previous paragraph we showed that the quantity associated with 
the shell-structured skeleton Q**^ should be constant for Fuclidean froths. In Fig. (4) {Q'"^^) / l-i2 is 
plotted (open symbols) in function of j for the different froths studied. One can note that in 
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all the systems analysed the quantity {Q'"^^) / l-i2 saturates at a value equal about to -2.0ib0.5 (S), 
-2.0ib0.2 (V), and -1.8ib0.2 (P). The saturation values should be similar for all three cases as 
{Q^^) 1 1^2 is a characteristic of the SSI skeleton. 

7 Conclusions 

We have studied the froth as organized in concentric layers of cells around a given central cell. 
Exact expressions for the number of cells in each layer, for the topological charge inside the shell 
cluster and for the average number of neighbors per cell in a given layer have been obtained. These 
topological properties of the shell-structure have been studied for the special class of SSI froths 
and for the general case where non-SSI defects are present. It turns out that the defects play a 
very important role in the organization of the froth structure. The defects enlarge the freedom for 
the construction of the cellular system around the germ cell by removing the constraints on the 
topological charge associated with Euclidean froths. With the relaxation of the constraints by the 
introduction of defects, we find a solution of space filling Euclidean froth which has only nearest 
neighbors correlations. In this solution, we have calculated approximate asymptotic expressions 
for /i'j, and which satisfy the Aboav relation and which are correlated only between first 
neighbors. Moreover, we have evaluated the average number of defects per layer. These expressions 
are free of adjustable parameters and describe well the behaviors of measurable properties of real 
froths and cellular patterns. 

Experimentally we find that soap froth in the steady state has an average number of cells per 
layer which grows linearly with the topological distance. The rate of growth is about 9. Slopes 
around 10 have been found for Voronoi froths. These slopes are considerably bigger than 2tt 
value suggested by simple geometrical consideration for the ratio between the perimeter of the 
shell-cluster and its radius. Moreover, we found that in soap and Voronoi froths the topological 
charge of the shell-cluster is always negative and decreases linearly with the size of the cluster. In 
particular, we observed that the slope is proportional to — /i2 and the coefficient of proportionality 
is smaller than 1 in soap froths and bigger than 1 in computer generated froths. The number of 
defects per layer have been also measured. Soap froths have in the asymptotic limit a percentage 
of defects around 10%, whereas bigger amounts have been found for computer generated froths. 

Our theory on the asymptotic behaviors is in good qualitative agreement with experiments. We 
correctly predict the linear growth of the number of cells per layer with a slope around 9. We 
demonstrate the linear decrement in the topological charge of the shell-cluster with a slope above 
the line defined by — /i2 in soap froths and with a slope below the line defined by — /i2 in Voronoi 
froths. We also predict percentages of defects per layer which are close to the experimental values 
and are smaller in soap froths and bigger in Voronoi froths. 

On the other hand, the quantitative agreement between the approximated predictions and the 
experimental data is not perfect. In the present paper we have obtained exact relations for the 
topological properties of the shell structure, but the predictions have been formulated under strong 
assumptions which simplify the exact results into expressions with no adjustable parameters. The 
partial disagreement between the approximated predictions and the experimental data might 
indicate that the assumption utilized are too strong or incorrect. Therefore there is still room for 
improvement. 
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A Generalization of the Weaire's sum rule 

Consider an z-cell of the froth with rii sides and the sum over the number of sides of the set of 
cells at a topological distance j from the z-cell. Such a sum is equal to mj[ni)Kj[ni) where mj[ni) 
denotes the average number of sides per cell in the layer distant j from the z-cell and Kj[ni) 
denotes the number of cells of this layer. Let us now sum this quantity over all the cells of the 
system: Y.- mj[ni)Kj[ni). In this sum the number of edges of each z-cell of the froth is counted a 
number of time equal to the number of cells at distance j from this cell (e.g. the number of sides 
of the generic z-cell is counted Kj[ni) times and contribute to the sum as niKj[ni). It follows that 
we have the identity 

= J2^^I''^ji^^) ■ (45) 

i i 

We can express this identity in term of the averages 

{m,K,) = (uK,) (46) 

where ((...)) indicates the average over the cell-sides distribution: ((...)) = J2n p{''^){---) with p[n) 
probability of an n-sided cell. 

Equation (46) is the generalization to the layer j of the Weaire's sum rule which is valid for the 
first layer (i.e {niin) = {n^))- Eor an arbitrary large system (arbitrary small boundary effects) 
relation (46) is exact and is valid also in the presence of non shell-reducible topological inclusions. 
One should note that, in real finite systems the effect of the boundary could be dramatically 
important. 



B Correlation and fluctuations 

A froth is uncorrelated after a given topological distance !^ if and only if for two cells respectively 
distant j > the probability Cj{n^m) to have one with n sides and the other with m sides 
factories Cj{n^m) = Sj[n)sj[m). It can be easily proved that in a froth where the cells are 
uncorrelated after the distance the average number of sides per cell (m"") in a layer j > ^ must 
be independent of the number of sides of the central cell. This is the physical consequence of the 
absence of correlations between the central cell and the cells in the layer j. This therefore means 
{mJ'^Kj) = m""'(/i'j). By using relation (46), we have 
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The topological charge inside layer j is 

- = (6 - m,)K, = ((6 - n)K,) + - r,_i , (48) 

The second equality comes from eq.(A2). For j > we put eq.(Bl) into (B2) to get for the 
uncorrelated case the relation 

r, = r,_i + ((6 - n)K,) - 1^ (for j > 0- (49) 

This is a recursive equation which gives the fluctuations Fj. Therefore, in uncorrelated froths 
the topological charge fluctuation are determined in term of the other statistical properties of the 
shell system. 

Consider for example a froth which is completely uncorrelated (topological gas, = 0). From 
eq.(49) one gets 

where we used the identities Ki = n and Fq = 6 — n. This relation gives the Aboav-Weaire's law 
(Fi = (1 — a)(6 — n), see in the text) with coefficient a = — /i2/6. 

C Proof for the correlation theorem of Euclidean SSI froth 



To prove this theorem, we first note that eq.(5) implies 

m = (^j-'E{Q^) ■ (51) 

8 = 

In order to fill the two-dimensional Euclidean space, there must be a minimum u such that 
(Kj) oc j for j > u. This implies that the cellular system must constrain the average topological 
charge inside a shell to be independent on the shell-size j. Such a constraint forces the average 
charge inside the layer j to be equal to zero for J > > 1, (see eq.l3) 

(g,) - = ((6 - m,)K,) = - X:((6 - n)F,) = - '-£{{6 - n)F,) = . (52) 

8 = 8 = 1 

Consequently, one gets the following two conditions on the fiuctuations F^ 

^((6-n)F,) = -/i2 forz/>l , (53) 

8 = 1 

and 

((6 - n)Fj) = for J > 1/ . (54) 

In appendix B we show that in a system where the topological correlations vanish after the 
layer, the fiuctuations must satisfy the following relation for j > ^ (eq.(B3)), 

((6-n)F,) = ((6-n)F,_i) + ii^^|^ . (55) 
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If one supposes f > ^ and if we set j = u m eq.(55), equations (52) and (54) yield ((6 — n)rj,_i) = 
which is in contradiction with the definition of u (for /i2 7^ 0). That implies u < ^. An immediate 
consequence is that in a random SSI froth where the Aboav-Weaire's law is satisfied with an 
arbitrary a (i.e. Fi = (1 — a)(6 — n), for some general a), the cells must be correlated at least 
between third neighbours. Indeed, u = 1 m eq.(54) implies (1 — a)/i2 = 0, or a = 1 if /i2 7^ 0, 
whereas u = 2 m eq.(53) implies (1 — a)/i2 = — /U2, or a = 2 if /i2 7^ 0. Both cases restrict a to 
have special values for /i2 7^ 0. It follows that if a is to remain a free parameter of the froth, 
u > 3. Since we have shown that ^ > therefore 1^ > 3. 
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aste et al Fig. 1 




Figure 1: Shell structure and defects in trivalent two-dimensional froth, the defects are shaded and 
the number denotes the topological distance from the center cell labeled (9, which is a deformed 
heptagon. 
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Aste et al Fig. 2 




Figure 2: Total number of cells per layer Kj in the j shell vs shell number j for soap (filled 
circle), random Voronoi construction (filled diamond) and Voronoi construction from perturbed 
triangular lattice (filled square). Total number of defects K'^ in the j-th shell vs j for soap (open 
circle), random Voronoi construction (open diamond) and Voronoi construction from perturbed 
triangular lattice (open square). 
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Aste et al Fig. 3 
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Figure 3: Defect concentration Aj vs shell number j for soap (open circle), random Voronoi 
construction (open diamond) and Voronoi construction from perturbed triangular lattice (open 
square). 
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Aste et al Fig. 4 
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Figure 4: Normalized cluster topological charge Qj/ f^2 vs shell number j for soap (filled circle), 
random Voronoi construction (filled diamond) and Voronoi construction from perturbed triangular 
lattice (filled square). The straight line indicate the slope -1. Soap and computer generated 
froths are in the opposite sides of this line. Normalized cluster topological charge of the skeleton, 
Q'j^ / 1-^2 vs shell number J for soap (open circle), random Voronoi construction (open diamond) and 
Voronoi construction from perturbed triangular lattice (open square). Note that the normalized 
topological charge for the skeleton saturated at about the same value for all cases. 



23 



